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ABSTRACT 

We have combined the well-studied CLASXS Chandra survey in Lockman 
with the Spitzer SWIRE survey data to study the X-ray-infrared connec- 
tion for AGN. The sample consists of 401 X-ray-sources, of which 306 are 
detected by Spitzer, and a further 257 AGN candidates detected through 
their dust torus, but not by Chandra. We have used spectroscopic redshifts 
and classifications from the literature, where available, and photometric 
redshifts for the remainder. For X-ray sources the X-ray hardness ratio 
has been modelled in terms of a power-law (F = 1.9) with absorption 
N(H). The optical and infrared data have been modelled in terms of our 
well-established optical galaxy and QSO templates, and infrared templates 
based on radiative transfer models. This type of analysis gives better in- 
sight into the infrared SEDs, and a better separation of the contribution 
of starbursts and AGN dust tori, than a simple comparison of 24 /im to 
optical or X-ray fluxes. We also believe this gives more insight than using 
a library of fixed UV- infrared templates. 

Our estimate of the N(H) distribution is consistent with other stud- 
ies, but we do find a higher proportion of low absorption objects at z 
< 0.5 than at z > 0.5. While we find only one X-ray AGN with N(H) 
> 10^"* cm~^, we argue that 10 objects with torus luminosity apparently 
exceeding the bolometric X-ray to 3/j,m luminosity are strong candidates 
for being heavily absorbed in X-rays. We also estimate that at least half of 
the infrared-detected AGN dust tori which are undetected in X-rays are 
likely to be Compton thick. Our estimate of the total number of Compton- 
thick objects in the 0.4 sq deg area is > 130, corresponding to > 20% of 
the combined SWIRE-CLASXS sample (and with an upper hmit of 39 %). 

We find no evidence for AGN with no dust tori, and none with a 
covering factor < 1% but there are clear examples of AGN with covering 
factors of only a few percent and these, though rare, do not fit easily with 
a unified picture for AGN. The range of dust covering factors is 1-100 
%, with a mean of 40%, ie a Type 2 fraction of 40%. Measured by the 
ratio of dust torus luminosity to X-ray or (for Type 1 objects) optical 
luminosity, the covering factor appears to decrease towards intermediate 
AGN luminosity, in contradiction to estimates based on ratios of narrow- 
line and broad-line spectra, but may increase again at low AGN luminosity. 

We find 7-10 candidate X-ray starbursts in the SWIRE-CLASXS sam- 
ple, with X-ray luminosities ranging up to Lxh = 10*^ erg s^^. This is a 
considerable extension of the luminosity range of X-ray starbursts previ- 
ously reported, but is consistent with the an extrapolation of the X-ray- 
infrared relation for starbursts into the realm of hyperluminous infrared 
galaxies. 

Key words: infrared: galaxies - galaxies: evolution - star:formation - 
galaxies: active - galaxies: starburst - X-ray:galaxies 
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1 INTRODUCTION 

The idea that AGN are surrounded by dust which 
absorbs their visible and ultraviolet light and reemits 
it at mid infrared wavelengths was first put forward 
by Rowan- Robinson (1976). He also suggested that 
extinction by this dust was responsible for the dis- 
tinction between Type 1 and 2 AGN, an idea which 
has developed into the unified model for AGN (An- 
tonucci 1993, Krolik 1999). IRAS showed that a mid- 
infrared excess is a common feature of AGN (Miley et 
al 1984) and this was interpreted in terms of a dust 
torus around the AGN nucleus (Rowan-Robinson and 
Crawford 1989, Pier and Krolik 1992, Granato and 
Danesc 1994, Efstathiou and Rowan- Robinson 1995). 
Rowan-Robinson (1995) proposed that the 'torus' 
might be an axisymmetric distribution of discrete, op- 
tically thick, clouds and gave a 1-dimensional treat- 
ment of the radiative transfer for such an ensemble 
of clouds. Recently 3-dirncnsional radiative transfer 
codes have been applied to such cloud ensembles by 
Hoenig (2006) and Nenkova et al (2008a,b). 

In this paper we explore the connection between 
X-ray emission from AGN and the mid-infrared emis- 
sion from the dust torus. We take advantage of the 
fact that the wide-area CLASXS Chandra survey 
(Yang et al 2004, Steffen et al 2004) lies within the 
SWIRE-Lockman area survey by Spitzor. This allows 
us to construct complete samples of AGN selected ei- 
ther in X-rays or at mid-infrared wavelengths, and 
to explore both the common sources and those un- 
detected in either waveband. The X-ray and optical 
properties of the CLASXS survey have been analyzed 
by Barger et al (2005). The X-ray mid- infrared con- 
nection has been explored in the SWIRE/Chandra 
survey of a 0.6 sq dog northern Lockman field, at 
similar X-ray sensitivity to CLASX, by PoUctta et al 
(2006), and in a shallower survey of 2.6 sq deg in the 
XMM MDS by Tajer et al (2007). Mid-infrared and 
X-ray spectral energy distributions (SEDs) of AGN in 
the XMM MDS have been analyzed by Franceschini et 
al (2005) and PoUetta et al (2007). In particular Pol- 
letta et al (2006) have discovered several Compton- 
thick AGN with prominent dust tori but very weaJi 
optical emission. Compton-thick quasars have been 
studied recently by Daddi et al (2007), Fiorc et al 
(2008, 2009) and Alexander et al (2008) (see also Pol- 
letta et al 2006 for references to earlier work). 

Our study is well-suited to address the issues of 
the prevalence of Compton thick AGN and whether 
there exists a major population of AGN which are 
absorbed even at 10 keV but might still contribute 
a significant fraction of the hard X-ray and infrared 
backgrounds. This is of interest because about half of 
the X-ray background is unresolved at 6 keV (Wors- 
ley et al 2005). Models have been proposed to explain 
the observed background which involve strong depen- 
dence of the obscured AGN fraction on X-ray lumi- 
nosity and redshift (eg La Franca et al 2005, Gilli et 



al 2007) or where the obscured fraction is indepen- 
dent of luminosity and redshift (eg Treister and Urry 
2005). Both types of model assume a high fraction 
of obscured AGN. We also explore whether there are 
AGN with very weak or absent dust tori and the role 
of X-ray starbursts in the sample. 



2 THE CLASXS CHANDRA SURVEY 

The CLASXS survey is described by Yang et al 
(2004). It covers an area of 0.4 sq deg in the Lock- 
man hole to a Chandra on-a:>ds X-ray sensitivity of 
5x10"^^ erg s"^ cm'^ at 0.4-2 keV and 3x10"^'^ erg 
s"^ cnr^ at 2-8 kcV . Optical photometry at BVRIz', 
spectroscopic redshifts and a spectroscopic classifica- 
tion are described by Steffen et al (2004) . For sources 
in common with SWIRE we also have UgriZ photom- 
etry. We have used the combined optical photometry 
to derive photometric redshifts, using the methodol- 
ogy of Rowan- Robinson et al (2008) , for sources with- 
out spectroscopic redshifts. Because we can be reason- 
ably certain that most of these sources are AGN, we 
have dropped the prior used in Rowan- Robinson et 
al (2008) that QSOs should be classified as optically 
starlike, and we have also allowed extinction values, 
Av up to 3.0 for QSOs. 

In the area in common with complete SWIRE 
coverage the CLASXS survey contains 424 hard X-ray 
sources of which 23 are stars identified spectroscopi- 
cally. 

2.1 Hydrogen column-density, N(H) 

We have followed Polletta et al (2006) in using the 
ratio of hard to soft X-ray count-rates to estimate 
N(H). We used PIMMS3 to derive a look-up table 
for the ratio of ACIS-I (2-8 keV) hard to soft (0.4- 
2.0) band counts as a a function of N(H) and red- 
shift, assuming an underlying power-law contirmum 
with r = 1.9 (Nandra and Pounds 1994). Tozzi et al 
(2006) give a mean value for F of 1.75, while Tueller 
et al (2008) give 1.98. Where a source is undetected 
in the soft X-ray band (6 sources), N(II) has been 
calculated assuming a soft count-rate of 0.02 and is 
a lower limit. However we should bear in mind that 
Winter et al (2008) find that 60% of a sample of hard 
X-ray sources have complex spectra, generally due to 
an additional soft X-ray excess, for which a simple 
power-law-|-absorption model will underestimate the 
true absorption. On the other hand Tozzi et al (2006), 
studying a sample similar to ours, find a smaller pro- 
portion (22/82 bright X-ray sources) showing either 
a scattered soft component or reflection-dominated 
spectrum. 

Figure 1 shows N(H) versus X-ray luminosity. 
The X-ray luminosity, Lxh, is the 4 keV rest-frame 
monochromatic luminosity, corrected for the assumed 
N(H). We have indicated the standard subdivisions 
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of N(H) into 'unabsorbed' (N(H) < lO^^cm-^), 'ob- 
scured' (lO^^cm^^ < N{H) < 10^^ cm"^), or 'Comp- 
ton thick' (N(H) > 10^* cm"^), and the 10*^ erg s"^ 
luminosity limit above which X-ray sources are almost 
certainly AGN. Figure 2L shows N(H) versus redsliift. 
There is a strong selection eflFect against detecting 
high values of N(H) at low redshift (see Fig 14 of Tozzi 
et al 2006). High values of N(H) (> lO^^-^cm"^) can 
only be detected for z>2, since at lower redshift the 
hard X-ray band would suffer serious absorption. The 
N(H)-z distribution seen in Fig 2L is similar to that 
found by Tozzi et al (2006). Fig 2R shows the same 
plot for X-ray sources with no SWIRE counterparts. 
The latter tend to be just a higher redsliift subsarnple 
of the X-ray population. 

Type 1 QSOs tend to have lower values of N(H), 
as expected in unified models, in which Type 1 QSOs 
are assumed to be viewed face-on. 

There are 88 sources for which the hardness ratio 
implies no absorption and these have been assigned 
a galactic absorption of 10^^'^® cm~^ (cf Yang et al 
2004). 

Figure 3 shows histograms of the N(H) distribu- 
tion for redshift bins between z = and 3. While 
the shifting of the high N(H) cutoff to higher values 
of N(H) as redshift increases is a selection effect as 
higher absorption drops sources out of the hard X-ray 
sample at lower z, the dearth of the low N(H) values 
at higher redshift seems to be a real effect. Because at 
higher redshift there is poorer resolution of low values 
of N(H), it is better to simply consider the total frac- 
tion of relatively unabsorbed sources. The percentage 
of low absorption (N(H) < lO^^crn^^) X-ray sources 
is 58, 37, 42, 32, 33, 39% for z = 0-0.5, 0.5-1, 1-1.5, 
1.5-2, 2-2.5, 2.5-3, respectively, so the increase at z 
< 0.5 seems real but there is no change in the frac- 
tion between z = 0.5 and 3. Tueller et al (2008) find 
for their low redshift Swift BAT 14-195 keV sample a 
flat distribution of logio N(H) from 20.5 to 23.5. Like 
Tozzi et al (2006), we seem to see fewer low N(H) sys- 
tems in the CLASXS Chandra sample than Tueller et 
al. 

2.2 X-ray bolometric correction 

Figure 4L shows Lopt versus Lxh for all the CLASXS 
sources. Fig 4R shows the corresponding plot just for 
objects with Type 1 QSO optical templates. For Type 
1 QSOs we can make a reasonably accurate estimate 
of the bolometric (X-ray to near ir) luminosity, since 
the SED peaks at UV wavelengths. We assume a bolo- 
metric correction to the 0.1-3 fmi luminosity, Lopt, of 
2.0 (Rowan- Robinson et al 2008). Where we have only 
X-ray luminosity as an indicator of the AGN luminos- 
ity, the issue is more problematic. In Fig 4R we see 
a wide range of ratios of optical to X-ray luminosity 
for Type 1 AGN, which translates to a wide range of 
bolometric corrections in the X-ray bands. Vasudevan 
and Fabian (2007) find bolometric corrections at 2-10 



keV in the range 4-100. Dependence of the bolomet- 
ric correction on luminosity have been investigated by 
Marconi et al (2004), Steffen et al (2006), Shemmer et 
al (2008), and on the Eddington fraction, L/LEdd, by 
Vasudevan and Fabian (2007). We find for the present 
sample (for optical QSOs with logioLxh > 43.0) a 
mean value of 27 for the bolometric correction in the 
hard X-ray (4 keV) band and have used this value to 
estimate the bolometric (X-ray to near ir) luminosity, 
Lxh,c, where necessary, but the large uncertainty in 
such estimates needs to be kept in mind. We have not 
attempted to correct for a dependence of bolometric 
correction on luminosity. 

We have investigated the outliers in Fig 4R. In 
Fig 5L we replot Fig 4R as the ratio Lopt,c/ Lxh,c ver- 
sus the ratio of hard to soft X-ray counts. QSOs with 
high values of Lopt/Lxh seem to be cases where the 
optical continuum may be contaminated by light from 
the parent galaxy, since they are flagged as optically 
extended. The QSO fits involve quite high values of 
Av (>1), and this suggests possible aliasing between 
highly reddened QSOs and galaxy SEDs in the optical 
template fitting. QSOs with high values of Lxh/Lopt 
are mostly high redshift QSOs, with modest redden- 
ing, but high X-ray absorption. They seem to be Type 
2 objects in which the optical light is scattered from 
the accretion disk into the line of sight. In Fig 5R we 
exclude objects flagged as optically extended or hav- 
ing N(H)> lCp"^an~^. There is now a much tighter 
correlation of optical and X-ray luminosities, with the 
corrected ratio Lopt,c/Lxh,c lying between 0.25 and 
10, still implying quite a wide range of hard X-ray 
bolometric corrections. 

Table 1(a) summarizes properties of the CLASXS 
X-ray sample in four bins of X-ray luminosity. 



3 SWIRE LOCKMAN SURVEY 

The Spitzer SWIRE survey is described by Lonsdale 
et al (2003, 2004). Source associations, spectral en- 
ergy distributions, photometric redshifts and infrared 
galaxy populations are discussed by Rowan- Robinson 
et al (2005, 2008). 

There are 8563 SWIRE sources from the SWIRE 
photometric catalogue (Rowan-Robinson et al 2008) 
in the CLASXS-SWIRE area. Of these, 233 are as- 
sociated with the 401 extragalactic CLASXS X-ray 
sources, using a 1.5" search radius between the respec- 
tive optical positions. To see if we can flnd SWIRE 
associations for the remaining X-ray sources we have 
examined the SWIRE sources that do not make it 
into the SWIRE photometric catalogue, either be- 
cause they are optically blank or do not meet the 
requirement of a total of two photometric bands out 
of UgriZ, 3.6, 4.5 fxm, of which at least one must be 
gri. A total of 2192 further SWIRE sources, satisfying 
either S(3.6) > 10 /uJy and S(4.5) > 8 /iJy, or S(24) 
> 100 /uJy, fall in the CLASXS-SWIRE area and 73 
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Figure 1. Hydrogen column-density, N(H), versus hard X-ray luminosity. The colour and symbol coding in this and 
subsequent figures are as follows: optical types (open symbols): black circle: optical QSO, green triangle: broad line, red 
square: galaxy, blue square: Seyfert 2, magenta square: no spectral classification, infrared template types (filled circles): 
black: cirrus, red: M82 starburst, green: Arp220 starburst, blue: AGN dust torus. The large filled circles are the two 
Compton-thick QSOs studied by PoUetta et al (2006). 



of these gave X-ray associations. So a total of 306 out 
of 401 extragalactic X-ray sources in CLASXS have 
SWIRE associations (76%). Although Spitzer detects 
a much larger population of galaxies than Chandra, 
the SWIRE catalogue still contains only 75% of the 
X-ray sources. The 96 non-SWIRE X-ray sources are 
shown in Fig 2R. They are typical of higher redshift 
X-ray sources and the non-detection by SWIRE does 
not imply unusual infrared properties. 

The number of SWIRE sources which are either 
Type 1 QSOs or have significant AGN dust tori com- 
ponents (contributing > 20% of the 8 or emis- 
sion) consists of 156 X-ray sources and 275 non-X- 
ray sources. So only 36% of the galaxies recognized 
as AGN by SWIRE, either through a QSO template 



in the optical or through the presence of an AGN 
dust torus, are detected in CLASXS. Thus both sur- 
veys are needed to fully characterize the AGN pop- 
ulation. SWIRE tends to miss some of the weaker 
dust tori, whereas CLASXSmisses X-ray absorbed or 
lower X-ray luminosity sources which SWIRE can de- 
tect through their dust tori. Donley et al (2005) first 
demonstrated the existence of AGN detected in the 
mid infrared by Spitzer through their dust tori, but 
undetected in deep Chandra integrations. Comparing 
a deep Chandra integration in the Groth strip with 
deep Spitzer data, Barmby et al (2006) found that 
90% of their 150 X-ray sources were detected by IRAC 
and two thirds by MIPS. 

If we regard 24 /im detection as essential for char- 
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Figure 2. LH: Hydrogen column-density, N(H), versus redshift. RH: Hydrogen column-density, N{H), versus hard X-ray 
luminosity, for X-ray sources with no SWIRE counterparts. 



acterization of dust emission, then 200 of the 401 ex- 
tragalactic CLASXS X-ray sources (50%) are SWIRE 
24 fim sources. Figure 6 shows the SWIRE 24 /im flux 
versus the hard X-ray flux for these sources. There 
is no detailed correlation and a very wide range of 
(S24/Sxh) ratio from ~0.004 to 1.0 (mJy/lO"^"* erg 
s~^ cm"'^) . The lowest values are candidates for very 
weak, low covering factor, dust tori, and high val- 
ues can indicate dust tori with very high covering 
factors, sources with X-ray absorption in the hard 
band which may be Compton thick sources, or possi- 
bly cases where X-ray emission is due to a starburst 
rather than an AGN. All these possibilities are dis- 
cussed below. 



3.1 Dust covering factor and Type 2 fraction 

Our ability to detect AGN dust tori at mid infrared 
wavelengths gives an excellent insight into the median 
dust covering factor for AGN. We fit infrared tem- 
plates only if there are at least two infrared bands 
with infrared excess relative to the starlight or QSO 
flt to optical data, and require that one of these bands 
must be 8 or 24 jj,m. Because AGN dust tori emit 
down to rest-frame wavelengths of 1 pm, in many 
cases dust tori are detected at all Spitzer wavelengths 
from 3.6-24 /xm (see Figs 11 and 12 below). Properties 
of infrared detected AGN dust tori are summarized 
in Table 1(b), binned by Ltor, and of the combined 
SWIRE-CLASXS AGN sample in Table 1(c). 

Figure 7L shows Ltor versus Lxh,c- Objects to 
the left of the locus of equal luminosity are either suf- 
fering from X-ray absorption or have an X-ray bolo- 
metric correction greater than the assumed value of 
27. We define the covering factor by the dust tori to 
be the ratio of the luminosity in the dust torus to the 



total bolometric luminosity of the AGN. From Fig 7L, 
the covering factor appears to range from 0.01 to 1. 

We define Lth, as the total (X-ray to 3 /im) bolo- 
metric luminosity of the AGN, excluding starburst or 
dust torus emission. For Type 1 QSOs, we assume a 
UV bolometric correction to the 0.1-3 /im luminosity, 
Lopt, of 2 (Rowan- Robinson et al 2008). Otherwise we 
estimate Lth from the X-ray luminosity, Lxh,c, using 
the bolometric correction of 27 derived for this sam- 
ple in section 2, bearing in mind that this will only be 
accurate on average. In the light of the discussion of 
Fig 5 above, we estimate Lth from X-ray luminosity 
even for Type 1 QSOs (i) if N(H) > 10^^ cm-^ (ii) if 
they are flagged as optically extended. 

Figure 8 shows Ltor /Lth versus Lx- We have in- 
dicated by vertical lines the luminosities below which 
X-ray sources may be due to starbursts, and above 
which they are almost certainly QSOs. The horizontal 
line at Ltor/ Lth ~ 1 corresponds to a dust torus cov- 
ering factor of 100%. Objects above this line are can- 
didates for being Compton thick AGN, though they 
could simply be cases of a higher than average X-ray 
bolometric correction if they are not Type 1 objects. 

The horizontal broken line corresponds to the 
mean dust torus covering factor of 40 % derived 
by Rowan- Robinson et al (2008) for a large sam- 
ple of SWIRE QSOs with Ltor > lO^'^ig, from < 
Ltor /Lopt >■ We can also estimate the mean covering 
factor for the present sample, from < Ltor/Lbh >, us- 
ing luminous X-ray QSOs {logwLxh > 44.5). We find 
a slightly lower mean covering factor of 35%, but this 
is based on only 9 QSOs and is not inconsistent with 
the 40 % derived by Rowan-Robinson et al (2008) us- 
ing 796 QSOs. The fraction of QSOs, as determined 
from optical template fitting, at logioLxh > 44.5 is 
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Figure 3. Histogram of the N(H) distribution for different redshift bins. Objects for which the inferred N(H) < 10^" cm ^ 

have all been assigned a value 10^^-^® 



25/46 (54%) again consistent with a 40 % covering 
factor. 

Objects with very low values of Ltor/Lhh are can- 
didate for having very weak or non-existent tori and 
are discussed in section 5 below. 

For sources with no X-ray detection we show 
Ltor/Lopt versus Ltor in Fig 9. We use Ltor as the 
X-axis because in all cases this is an, albeit indi- 
rect, measure of the AGN luminosity. We have plot- 
ted this for SWIRE AGN dust tori with no hard X- 
ray detection (Fig 9L) , and the hard X-ray AGN 
(Fig 9R), separately. For X-ray sources the distribu- 
tion is broadly equivalent to Fig 8, except that lines 
of constant Lopt now have slope (we have indi- 
cated with a broken line the approximate boundary 
between galaxies and QSOs). AGN which have opti- 
cal QSO templates and/or broad lines show a trend 
in which the covering factor {Ltor/Lopt) decreases as 



Ltor decreases. Such a trend can also be seen in Fig 

8. For AGN dust tori associated with galaxies (the 
overwhelming majority of the non-X-ray AGN), the 
sources may either be Type 2 objects, in which the 
QSO is shrouded by dust, or weaker Seyferts where 
we do not see the optical AGN against the light of 
the host galaxy. Ltor /Lopt will underestimate the true 
covering factor for Type 1 Seyferts. However for ob- 
jects with Ltor > 2Lopt we can be confident they are 
Type 2 objects, since even for an almost 100% cov- 
ering factor the QSO would outshine the galaxy if 
viewed pole-on. LtorfiLopt clearly overestimates the 
covering factor for these objects. Despite these limi- 
tations it is of interest to see how the median value of 
Ltor/^Lopt for the combined infrared and X-ray sam- 
ple varies with torus luminosity and we have indicated 
this iu the last two lines of Table 1 (c). We find that 
Ltor/2Lopt decrease reduced, from 0.4 in 
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Figure 4. Optical luminosity, Lopt, versus X-ray luminosity, • LH plot: all sources, RH plot: optical QSOs only. The 
solid line corresponds to equal bolometric luminosities assuming an optical-UV bolometric correction of 2.0 and a hard 
X-ray correction of 27. Crosses denote Compton-thick sources {N(H) > 10^'*cm~^). 
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Figure 5. Ratio of optical luminosity, Lopt, to X-ray luminosity, Lxi versus ratio of hard to soft X-ray counts. LH plot: all 
sources assigned QSO templates, RH plot: Same, excluding QSOs flagged as optically extended or with N(H)> 10^'' cm"^. 
The solid lines correspond to apx = 1-3, 1.7. The broken line corresponds to an unabsorbed F =1.9 power-law. 



the highest luminosity bin, consistent with the result 
of Rowan- Robinson et al (2008), to 0.08 in the lowest 
luminosity bin (this trend is indicated by the dotted 
line in Fig 9). However the lowest luminosity dust tori 
are probably simply Seyferts, AGN with optical lumi- 
nosities much less than that of their parent galaxy, so 
Ltor/2Lopt severely underestimates the covering dust 
factor. In Table 1(a) we have shown the median value 
of lg{Ltor/Li,h) for the X-ray sample. Our data are 
consistent with a model in which the median cover- 
ing factor at first decreases from high to moderate 
luminosity, then increases towards low luminosities. 



This would run counter to several claims in the liter- 
ature, which suggest that the covering factor increases 
monotonically towards lower luminosities. 

Hasinger (2008) showed that the ratio of narrow- 
line X-ray AGN to broad-line AGN increases with 
decreasing hard X-ray luminosity, using nine X-ray 
surveys of which CLASXS is one. A similar trend is 
found also by e.g. Ueda et al. (2003), La Franca et 
al. (2005), Gilh et al. (2007), and Fiore et al. (2008). 
However all these studies use the fraction of absorbed 
AGN (from X-rays) or broad line AGN (from opti- 
cal spectroscopy) which give information on the total 
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line-of-sight absorption, which may include galaxy- 
scale absorption (for Seyfert 1.8 andl.9, Maoilino and 
Rieke 1995) or dust-free ionized gas (Comastri et al 
2001, Nandra et al 2003). In Table 1 we show the 
corresponding numbers of narrow-line and broad-line 
AGN in our sample, and their ratio, which does in- 
deed show this effect. However there are some prob- 
lems with immediately interpreting this as an in- 
crease of covering factor with declining luminosity. 
The presence of broad lines, or of an optical QSO 
continuum, is not a reliable indicator that an ob- 
ject is face-on. Several of the objects that we believe 
have to be Type 2, Compton-thick AGN, on the ba- 
sis of Ltor/Lth or Ltor/Lopt, have broad-line spec- 
tra and/or an optical QSO continuum (see section 4 
and Table 2). PoUetta et al (2006) 's object SWIRE- 
J104409.95+585224.8 has Ltor ~ WOLopt, yet its op- 
tical spectrum is that of an unreddened broad-line 
QSO, clearly because it is seen in scattered light in 



the optical. A second problem with arguments based 
on broad/narrow line ratios is that they are based on 
a small fraction of the X-ray AGN sample (22% of the 
AGN with Lxh < W^^'^ erg s'^ in Table 1(a)). If we 
look more closely at the lowest luminosity bin, where 
5 out of 7 AGN are narrow-line, it turns out that 4 of 
these are low absorption AGN (the 5th has Zo(7ioN(H) 
— 22.38) and for 3 of them their optical SEDs are fit- 
ted with QSO templates. So these 5 objects do not on 
their own constitute a strong case for a high covering 
factor at low luminosities. 

To investigate this further we have divided all 
CLASXS objects spectroscopically classified as broad- 
line or narrow-line into unabsorbed [logioN [H) < 
22), absorbed (22 < logioN{H) < 23) and strongly 
absorbed {logioN{H) > 23). In the unified model for 
AGN we might have expected that all broad-line ob- 
jects would be unabsorbed and that all narrow-line 
objects would show at least some absorption. The ac- 
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Figure 7. LH: AGN dust torus luminosity, Ltor, versus bolometric X-ray luminosity, Lxh c- RH: Starburst luminosity, 
Lsb, versus bolometric X-ray luminosity, Lxh,c- The broken lines denote the Ranalli et al (2003) relation for X-ray emission 
from star bursts. 



tual numbers of broad-line objects in the 3 classes 
are 53, 31 and 5 respectively, and the numbers of 
narrow-line objects are 15, 15 and 7 respectively. So 
the broad/narrovvf-line distinction is not a very good 
indicator of X-ray absorption. The broad-line objects 
with strong X-ray absorption presumably may imply 
that the broad lines are scattered round the X-ray ab- 
sorbing cloud or clouds. Alternatively there may be 
a link with BAL QSOs which are generally not very 
strongly absorbed in the optical but are characterized 
by strong X-ray absorption (Gallagher et al 2006, Fan 
et al 2009), due to a dust-free outflowing medium. To 
explain the narrow-line unabsorbed objects we must 
presumably invoke excess soft X-ray emission mask- 
ing the absorption. 

Maiolino et al (2007) have presented evidence of 
increasing Type 2 fraction as luminosity is reduced, 
based on the ratio of the 6.7 /im to the 5100 A contin- 
uum, which should be equivalent to our Ltor / Lopt ra- 
tio. However this is not supported by the much larger 
and perhaps more representative sample of our Fig 9. 
Treister et al (2008) have also argued for an increas- 
ing dust covering factor towards lower luminosities 
by considering the ratio of 24 fim luminosity to bolo- 
metric luminosity. We believe it is crucial to separate 
the torus and starburst contributions to the 24 /im 
emission. 

Some authors (eg La Franca et al (2005), 
Hasinger 2008, Treister et al 2008) have suggested 
that the fraction of obscured AGN increases signifi- 
cantly with redshift. Because higher luminosity AGN 
tend to be at higher redshifts in flux-limited samples 
like ours, this could account for part of the increase 
in covering factor seen for the highest luminosity bin 
in Table 1. Our sample is probably not large enough 
to disentangle luminosity and redshift dependence. 



3.2 Fraction of Compton-thick AGN 

Figure 10 shows Ltor versus redshift (LH) and Lxh 
versus redshift (RH), illustrating the much greater 
sensitivity of SWIRE to lower-luminosity, lower-z, 
dust tori, while Chandra sees deeper in the high- 
luminosity sources. It is apparent that many of the 
lower luminosity dust tori would by easily detectable 
in X-rays if they had similar properties to the X-ray 
detected objects. To quantify this we have estimated 
the unabsorbed X-ray luminosity each of these sources 
would have with a hard X-ray flux at the limit of the 
CLASXS survey {3.10~^^ergs~^cm~^), and then used 
the ansatz 

Lth ~ 27Lxh ~ 'iLopt ~ 2.51/tor (1) 
appropriate to a dust torus with the average op- 
tical and X-ray bolometric corrections, and average 
dust covering factor of 40%, to estimate how much 
X-ray absorption would be needed to push the X-ray 
flux down to this limit. 116 out 246 infrared AGN dust 
tori not detected in X-rays were found to need a min- 
imum hydrogen column N(H) of ~ 10^** cm~^, cor- 
responding to Compton-thick objects. Table 1 sum- 
marizes some properties of the X-ray and non-X-ray 
sources in the SWIRE-CLASXS sample, as a function 
of luminosity class. Our conclusion is that > 20% of 
the sample is Compton thick, and that > 52% have 
N(H) > 10^^ cm~'^, but we are not able to demonstrate 
any clear dependence of these percentages on X-ray 
or torus luminosity. ~ 25% of the sample have N(H) 
< lO^^cm"^ and are unabsorbed. Again there is little 
dependence on AGN luminosity except perhaps that 
the the proportion of unabsorbed sources is higher in 
the lowest bin of X-ray or torus luminosity. Our re- 
sults therefore favour the X-ray background model of 
Treister and Urry (2005), who assume a 25 % frac- 
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Figure 8. Ratio of torus luminosity to bolometric optical-X-ray luminosity versus X-ray luminosity. Crosses denote upper 
limits for AGN with no dust torus detected by SWIRE. 



tion of unabsorbed sources independently of luminos- 
ity and of redshift. However we could not rule out a 
model in which the Compton-thick fraction went from 
20% at Lxh = lO*^-^ erg s'^ to 50% at Lxh = iC^'^ 
erg s~^, and such a behaviour might contribute to a 
viable model of the hard X-ray background. 



4 MODELLING THE SEDS OF 
COMPTON THICK AGN 

The optical-infrared SED of the one Compton thick 
AGN identified from the ratio of hard to soft counts 
is shown in Fig IIL (object 1, Table 2a). It is a Type 
1 QSO in its optical SED, and has Ltor > Lopt, but it 
does not have the vey high Ltor /Lopt ratio of the two 
best studied Compton-thick objects of PoUetta et al 
(2006). There may be more Compton-thick objects in 
the CLASXS X-ray sample because of the possibility 



of excess soft X-ray emission from ionized gas (warm 
absorber model, Gierlinski and Done 2006), Compton 
scattered fight (refiected from accretion disk, Nandra 
and Pounds 1994, Maloney and Reynolds 2000, Ross 
and Fabian 2005, Nandra and Iwasawa 2007) or inho- 
mogenous covering (Winter et al 2008). The latter is 
to be expected in a multi-cloud torus model. 

Fig 11 also shows 9 candidate Compton thick X- 
ray objects from the condition Ltor/Lbh > 1- By defi- 
nition all have prominent dust tori. Because of uncer- 
tainties in the bolometric correction factors used to 
calculate Lth, we can not be certain that any of these 
objects require enhanced X-ray absorption. A further 
uncertainty could be introduced by the possibility of 
X-ray variability in the interval between the Chan- 
dra and Spitzer observations. The direct estimates of 
N(H) are in the range 10^^ - 4.10^^, and 3 are unab- 
sorbed, but these may be underestimates because of 
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Table 1. Properties of X-ray and infrared AGN in SWIRE-CLASXS sample, as a function of luminosity 



property 




]^q10.35 


< Lfc^ < 1011-35 


1011-35 < < 1012-35 


Lbh > 10^2-35 


total 
















(a) X-ray detected AGN 














Lxh (erg s-i) 


< 1042-5 


1042-5 . 


- 1043-5 


1043-5 _ 1044.5 


> 1044-5 




tot 3.1 numbGr 




117 




1 




401 


X-ray starbursts 


4 


5 




1 





10 


optical QSO template 


9 


21 




51 


25 


106 


IN in. 1 ^ j_u c//fc 


n 
u 


u 




u 


1 


1 


1023 <N(H)< 1024 





8 




61 


23 


92 


N(H) < 1022cm-2 


34 


63 




53 


9 


159 


Ltor > Li,ti 





4 




8 


1 


13 


no. narrow-line AGN 


5 


15 




14 


3 


37 


no. broad-line AGN 


2 


12 




57 


18 


89 






1 2 






17 
U. 1 / 


0.42 


median lg(Ltor / Li,h) 


-0.60 


-0.75 




-0.80 


-0.58 




implied covering factor 


25% 


18% 




16% 


26% 




(b) infrared detected AGN dust tori 














(no X) Ltor{LQ) 


< 109-95 


109-95 - 


_ 1010-95 


1010.96 _ 1011-96 


> 1011-95 




total number 


56 


175 




18 


8 


257 


optical QSO template 













7 


7 


N(H) > 1024cm-2 


>o 


>92 




16 


8 


>116 


N(H) > 1023 cm-2 


47 


174 




18 


8 


247 




1 


12 




3 


1 


17 


(c) combined sample 


96 


292 




216 


54 


658 


(excluding starbursts) 














no. optical QSO template 


>9 


>21 




>51 


32 


>114 


% QSO 


> 8% 


> 7% 




> 24% 


60% 




no. Compton-tliick 


>0 


>96 




24 


10 


>130 


% Compton-tliick 


>0% 


> 33% 




> 11% 


18% 


> 209 


no. N(H) > 1023cm-2 


47 


182 




79 


32 


247 


% N(H) > 1023 cm-2 


50 


63 




39 


58 


41 


median lg(Ltor/2Lopt) 


-1.1 


-0.55 




-0.7 


-0.4 




implied covering factor 


8% 


28% 




20% 


40% 





the cfTccts of Compton scattering. Objects 3, 5, 6, 10 
in Table 2b and Fig 11, which have Ltor » Lxh,c and 
optical seds fitted with galaxy templates, are candi- 
dates for being obscured AGN of the type discussed 
by Ueda et al (2008), where the torus covering factor 
is close to 100%. 

A further 18 SWIRE objects not detected in X 
rays have Ltor > Lopt,c, but 17 of these are galax- 
ies with much lower redshifts and torus luminosities 
than the objects plotted in Fig 11 and could be simply 
Type 2 QSOs. However the estimate of the absorp- 
tion required to explain their non-detection in X-rays 
(see section 3) shows that most are probably Comp- 
ton thick, and all have N(H) > 102^ * cm-2. rpj^^ -^g^j^ 
object has a QSO optical SED, with photometric red- 
shift 2.89, and is very similar to the X-ray objects in 
Fig 11. It has been included in Fig 11 (object 11). The 
detailed model of Fig 11 has Ltor < Lopt,c but the es- 
timate of the absorption required to explain its non- 
detection in X-rays shows that it is definitely Comp- 
ton thick. Properties of both types of Compton thick 
candidates are given in Table 2. Spectroscopic types 



in this and subsequent tables are: 4 = broad-line QSO, 
2,3 = narrow-line AGN, 1 = galaxy. Photometric red- 
shifts are shown bracketed. 



According to the AGN dust torus models of Efs- 
tathiou and Rowan-Robinson (1995), a highly edge-on 
dust torus would have its emission pealcing at a rest- 
frame wavelength ~ SO/um. We have not seen any such 
objects but do not have enough sensitivity at 70 and 
160 /im to readily detect such tori. 



There are 73 X-ray sources in the SWIRE)- 
CLASXS sample with z>2, and the number of 
Compton-thick candidates among these is one directly 
from N(H), 4 from Ltor > Lopt,c, with perhaps a 
further 10 objects missed from the sample due to 
strong X-ray absorption (section 3) , corresponding to 
a Compton-thick fraction 15/83 or 19% (see also Ta- 
ble 1). 
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Figure 9. Ratio of torus luminosity to optical luminosity versus torus luminosity, for: (L) infrared-selected AGN (SWIRE 
sources with optical QSO templates, or presence of AGN dust tori) and no hard X-ray detection, (R) Hard X-ray AGN 
sources with optical QSO templates, or presence of AGN dust tori . The large filled circles are the two Compton-thick 
QSOs studied by PoUetta et al (2006). The broken line indicates the optical luminosity above which objects are likely to 
be QSOs. The dotted line indicates the trend of the median Ltor/ Lopt for the whole SWIRE-CLASXS sample (last line of 
Table 1). 




Figure 10. LH: Torus luminosity versus redshift. Filled red circles: starburst ir SEDs, filled blue circles: dost torus 
dominated; green open triangle: broad-line; blue open square: narow-line; red open square: galaxy spectra; black open 
circle: optical QSO. RH: X-ray luminosity versus redshift. Sources with infrared dust tori: open green circle; and with 
optical QSO (filled blue circle); with E-Sab (filled green circle); with Sbc-sb (filled red circle). Other X-ray galaxies: filled 
black circle. 



5 WEAK TORI AGN 

Figure 12L shows the SEDs of X-ray AGN which are 
candidates for having very weak tori (properties given 
in Table 3). For five objects (2, 4, 5, 6, 10) the pho- 
tometric redshift solution selected a highly reddened 
QSO, which then gives values for Lth factors of 5-20 
times Lxh,c- However they are all flagged as optically 



extended objects and galaxy models seem much more 
plausible. The ratios Ltor/Lxh,c become 0.2-0.3. The 
remaining objects (1, 3, 7, 8, 9) appear to be Type 
2 objects, with Lxh,c > Lopt- These show a range of 
covering factors from 1-3%. 

There were a further 3 objects which our auto- 
matic SED-fitting code gave strongly reddened QSO 
fits in the optical, with Arp 220 or cirrus+starburst 
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Table 2. Compton thick candidates: (a) N{H) > lO^^cm-^, (b) Ltor > L^h 



object no. 


RA 


dec 


z 


spec, type 


■^xhc 


logioN{H) 




Ltor 


Lopt 


type 


(a) 
1 


158.55830 


57.77818 


(2.68) 




13.14 


24.04 




12.27 


12.17 


QSO 


(b) 






















2 


158.34232 


57.54998 


2.096 


4 


11.33 


19.57 


12.21 


11.37 


12.03 


sb 


3 


158.68903 


57.67629 


(1.87) 




10.68 


19.57 


12.23 


11.33 


11.61 


Sbc 


4 


158.33299 


57.80566 


1.928 


4 


12.05 


21.13 


12.48 


12.24 


12.12 


sb 


5 


158.78539 


57.96220 


1.261 


4 


11.84 


23.19 




12.45 


12.65 


sb 


6 


158.92015 


57.84380 


1.370 


4 


11.56 


22.62 


12.20 


12.14 


12.34 


sb 


7 


158.32928 


57.85759 


1.823 


4 


12.04 


22.65 


12.35 


12.14 


12.32 


sb 


8 


158.30424 


57.56830 


(2.565) 




12.24 


22.85 




12.33 


10.61 


Scd 


9 


158.41524 


57.63800 


(1.366) 




10.91 


22.67 


11.88 


11.16 


9.70 


Scd 


10 


158.72195 


57.56465 


1.270 




10.81 


19.57 




11.50 


9.95 


Scd 


11 


158.45297 


57.70647 


(2.63) 




12.15 


23.61 




12.02 


11.07 


QSO 


12 


158.66422 


57.59955 


(2.89) 








12.90 


12.92 


13.32 


QSO 



fits in tlie infrared, ie no dust torus. However more 
careful modeling of the SEDs showed that these were 
probably unreddened galaxy SEDs in the optical, with 
dust tori and starbursts in the infrared, ie they are in 
fact Type 2 objects. Thus there are no clear cases 
where the implied limit on the dust torus component 
would correspond to a covering factor < 1%. Our con- 
clusion is that the dust covering factor around lumi- 
nous AGN ranges from 1-100%, with a median value 
of 40%. However dust covering factors of just a few 
percent are not very different from naked quasars and 
these objects with very weak tori do not fit easily into 
a unification model in which AGN properties depend 
only on viewing angle. 



6 THE AGN-STARBURST CONNECTION 
AND X-RAY STARBURSTS 

X-ray emission from starbursts, due to associated X- 
ray binaries, supernova remnants and difiiuse hot gas, 
have been studied by Griffiths and Padovani (1990), 
David et al (1992), Nandra ct al (2002), and Ranalh 
et al (2003). Ranalli et al (2003) use 23 local starburst 
galaxies to derive a mean ratio of far infrared to X-ray 
luminosity for staxbursts. Namdra et al (2002) show 
that this ratio is very similar for z ~ 3 Lyman-break 
galaxies. 

Figure 7R shows Lsb versus Lx for our SWIREl- 
CLASXS sample. We have included A220 and cirrus 
dominated sources in this plot, with Lir replacing Lsb- 
The range of ratios derived by Ranalli et al (2003) for 
X-ray emission from starbursts is shown as a pair of 
broken lines. Five sources lie within this range, with a 
few others just below the lower edge. Fig 12R shows 
SEDs for sources with logioLir > logioLxh,c + 1-3 
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Figure 12. LH: SEDs of objects with very weaJj tori (Table 3). RH: SEDs of candidates for being X-ray starbursts (Table 
4). 



Table 3. X-ray sources with weak dust tori 



object no. 


RA 


dec 


z 


spec, type 


^xhc 


logioN{H) 






Lopt 


type 


Av 


1 


158.81332 


57.51605 


(0.959) 




12.39 


23.41 


11.75 


10.60 


11.44 


Sbc 


0.4 


2 


158.73761 


57.55341 


0.937 


3 


10.93 


19.57 


11.22 


10.47 


11.32 


Scd 


0.2 




















(11.95 


QSO 


1.4) 


3 


158.74062 


57.63224 


(0.950) 




11.50 


23.07 


11.61 


9.52 


11.14 


Scd 


1.4 


4 


158.66597 


57.63459 


0.503 


4 


9.74 


19.57 


10.56 


10.04 


10.74 


Scd 


0.2 




















(11.41 


QSO 


1.3) 


5 


158.45320 


57.69691 


(0.282) 




9.69 


21.28 


9.79 


8.74 


9.49 


Sbc 


0.6 




















(10.45 


QSO 


2.4) 


6 


158.66617 


57.73209 


0.452 


1 


11.37 


21.40 


11.28 


10.43 


11.13 


Sbc 


0.0 




















(12.00 


QSO 


1.6) 


7 


158.99966 


57.73656 


0.516 


2 


11.94 


22.81 




10.15 


10.85 


Sbc 


0.3 


8 


158.69876 


57.83265 


1.015 


3 


11.81 


22.76 


11.93 


10.03 


11.43 


Sab 




9 


158.88411 


57.94242 


0.947 


1 


11.87 


23.12 


11.12 


10.37 


10.72 


E 






















10.4 


QSOl 




10 


158.35919 


58.02242 


0.736 


4 


11.55 


21.09 


11.72 


10.75 


11.46 


Scd 






















(11.91 


QSOl 


0.7) 



which have 24 /im detection (properties given in Ta- 
ble 4). Apart from object 8, which shows some ev- 
idence for dust torus emission, all the other objects 
are good candidates to be X-ray starbursts. All ex- 
cept object 4 would be classified as unabsorbed X-ray 
sources. While 4 of the 10 objects have low X-ray lu- 
minosities {Lxh < W^^ergs~^), the rest are of higher 
luminosity. Three are hyperluminous infrared galax- 
ies and have correspondingly strong X-ray emission. 
It appears that X-ray starbursts can have X-ray lumi- 
nosities ranging up to lO^^ergs"^. Spectroscopically 
object 1 is classifed as broad-line and object 5 is clas- 
sified as narrow-line, but neither shows evidence for 
an optical QSO or for a dust torus. We conclude that 
at least 7 of these objects are X-ray starbursts, with 
objects 1, 5 and 8 showing some evidence for the pres- 
ence of an AGN. It would be valuable, especially for 



the starburst candidates with z>l, to obtain spectro- 
scopic redshifts for them. 



For the rest of the CLASXS sample the X-ray 
emission is presumably due to an AGN. The strong 
correlation in Fig 7R is therefore interesting, suggest- 
ing a common gas feeding mechanism for both star- 
bursts and for black hole accretion, presumably due to 
galaxy interactions and mergers. Since Lsb is propor- 
tional to the star formation rate, and for a constant 
Eddington factor, Lxh is proportional to the black 
hole accretion rate, this plot can be seen as the time 
derivative of the Magorrian relation between bulge 
stellar mass and black hole mass. An earlier version 
of this correlation, in the form Lsb versus Ltor, was 
given by Rowan- Robinson (2000) . 
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Table 4. X-ray starburst candidates 



object no. 


RA 


dec 


z 


spec, type 


^xhc 


1 


158.67293 


57.54453 


1.208 


4 


11.10 


2 


158.67242 


57.89325 


(1.570) 




11.99 


3 


158.49472 


57.72141 


0.264 




9.71 


4 


158.85077 


57.74316 


0.234 


1 


9.68 


5 


158.54678 


57.92445 


0.203 


2 


9.08 


6 


158.32266 


57.87665 


0.094 


1 


8.07 


7 


158.14431 


57.93956 


(1.432) 




10.94 


8 


158.68903 


57.67629 


(1.871) 




10.68 


9 


158.21892 


57.74105 


(1.559) 




10.74 


10 


158.71025 


57.68777 


(2.819) 




11.29 



logioNiH) 






Ltor ^opt 


type 


Av 


21.65 




12.68 (A220) 


11.12 


Sod 


0.0 


19.57 




13.34 (A220) 


11.60 


sb 


0.4 


21.61 


11.16 


10.81 


11.64 


Sod 


0.9 


22.27 




10.87 


10.89 


Sod 


0.9 


21.04 




10.96 


10.56 


Sod 


0.6 


19.57 


9.16 




9.40 


Scd 


0.0 


19.57 




13.36 (A220) 


10.91 


Scd 


0.8 


19.57 




12.23 


11.33 11.49 


Sbc 


0.0 


19.57 




12.30 


11.82 


sb 


0.3 


19.57 




13.17 


11.92 


Sbc 


0.4 








(11.5 


QSO 


0.7) 



7 DISCUSSION AND CONCLUSIONS 

As has been argued in several recent papers (Polletta 
et al 2006, Daddi et al 2007, Fiore et al 2008, 2009), 
combining data from X-ray surveys and from Spitzer 
gives us a better understanding of the statistics of the 
dust-covering factor around AGN. We have combined 
the well-studied CLASXS Chandra survey in Lock- 
man with the SWIRE survey data in the overlapping 
area of the sky, 0.4 sq deg. The sample consists of 401 
X-ray-sources, of which 306 axe detected by Spitzer, 
and a further 257 AGN detected through their dust 
torus, but not by Chandra. We have used the spec- 
troscopic redshifts and classifications of Steffen et al 
(2004), and other spectroscopy from the literature, 
where available, and photometric redshifts from the 
methodology of Rowan-Robinson et al (2008) for the 
remainder. For X-ray sources the X-ray hardness ra- 
tio has been modelled in terms of a power-law (F = 
1.9) with absorption N(H). The optical and infrared 
data have been modelled in terms of the galaxy and 
QSO templates, and infrared templates, of Rowan- 
Robinson et al (2008). The present analysis in terms 
of well-established infrared templates based on radia- 
tive transfer models gives better insight into the in- 
frared SEDs, and a better separation of the contribu- 
tion of starbursts and AGN dust tori, than a simple 
comparison of 24 /im to optical or X-ray fluxes. We 
also believe this gives more insight than using a li- 
brary of fixed UV-infrared templates. 

Our estimate of the N(H) distribution is consis- 
tent with other studies, but we do find a higher pro- 
portion of low absorption objects at z < 0.5 (58% 
of the X-ray sample) than at z > 0.5 (37%). X-ray 
selection effects make it difficult to get a good esti- 
mate of the total number of high extinction AGN, 
especially at z < 2. Almost half of the AGN dust 
tori which we detect with Spitzer but not with Chan- 
dra arc likely to be Compton thick objects and we 
estimate that at least 20% of our combined AGN 
sample are Compton-thick objects. If all the sources 
detected with Spitzer but not Chandra were Comp- 
ton thick objects, which we regard as highly unlik- 
ley, the overall percentage of Compton-thick objects 



would be increased to 39%. We could not however 
rule out a model in which the Compton-thick frac- 
tion went from 20% at Lxh = lO"*"*'^ erg s"^ to 50% 
at Lxh = lO''^'® erg s~^, and such a model might pro- 
vide a viable model of the X-ray background (Gilli et 
al 2007). 

A further uncertainty for the X-ray sample arises 
from the possible presence of enhanced soft X-ray 
emission (Nandra and Pounds 1994, Maloney and 
Reynolds 2000, Ross and Fabian 2005, Gierlinski and 
Done 2006, Nandra and Iwasawa 2007, Winter et al 
2008) , which can mask the presence of Compton-thick 
objects. We suggest that the objects in Table 2b fall 
into this category. 

7.1 Comparison with other work on 
Compton thick AGN 

Polletta et al (2006) identify 7 X-ray sources with 
N(H) > 10^* cm~^, estimate a total of 55 in 0.6 sq deg 
of Chandra/SWIRE area, only 20% of which are de- 
tected in X-rays. While we find only one X-ray AGN 
with N(H) > 10^* cm~^, our estimate of the total 
number of Compton-thick objects in our 0.4 sq deg 
area, to comparable Chandra depth, is > 130, corre- 
sponding to > 20% of the combined SWIRE-CLASXS 
sample. 

Daddi et al (2007) note that 20-30% of faint z= 
1.4-2.5 galaxies detected with MIPS at 24 jira have 
mid-infrared excess unlikely to be due to star for- 
mation, stacked X-ray spectra rising steeply at > 
10 KeV, suggesting they host Compton thick AGN, 
and have space-densities twice that of X-ray detected 
AGN. This is very consistent with our estimate of 
> 20% for a sample most of whose redshifts range 
from 0.4-3. 

Fiore et al (2008) studied a population of objects 
in Chandra DFS with high 24 ^m to optical ratios 
and estimated that 80% likely to be Compton thick, 
with the number at z = 1.2-2.6 similar to the number 
of unobscured and moderately obscured AGNs. Fiore 
et al (2009) estimate that 44% of X-ray sources with 
logioLxh =44-45, and 67% with logwLxh = 43.5- 
44, are Compton thick, from COSMOS MIPS survey. 
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These are much higher percentages than we find (see 
Table 1). 

Alexander et al (2008) note the rarity of reported 
Compton-thick objects at z = 2-2.5. They discuss a 
z = 2.211 Cornpton thick quasar for which they have 
optical and IRS spectroscopy and identifiy a further 6 
at z = 2-2.5 in the literature. The present study identi- 
fies a further 4 X-ray selected Compton-thick objects 
at z > 2 on the basis of photometric redshifts (Table 
2), and there arc a further 6 infrarcd-dctcctcd AGN at 
z > 2 which we estimate to be Compton-thick. Clearly 
it is desirable to obtain optical spectroscopy for these 
Compton-thick candidates. 

7.2 AGN with weak tori, X-ray starbursts 

We conclude that there is no evidence for AGN with 
no dust tori, and none with a covering factor < 1%. 
The range of dust covering factors is 1-100 %, with a 
mean of 40%, ie a Type 2 fraction of 40%. However 
dust covering factors of just a few percent are not 
very different from naked quasars and these objects 
with very weak tori do not fit easily into a unifica- 
tion model in which AGN properties depend only on 
viewing angle. 

Measured by the ratio of dust torus luminosity 
to X-ray or (for Type 1 objects) optical luminosity, 
the covering factor appears to decrease towards lower 
AGN luminosity, and certainly shows no evidence of 
an increase, in contradiction to estimates based on ra- 
tios of narrow-line and broad-line spectra. The lines- 
of-sight to X-ray emitting, line emitting and optical 
continuum emission are complex and often involve 
scattering processes. The material responsible for the 
soft X-ray absorption is not necessarily well correlated 
with the dust responsible for the dust torus emission. 

We find 7-10 X-ray starbursts in the SWIRE- 
CLASXS sample, with X-ray luminosities ranging up 
to Lxh = 10** erg s^^. This is a considerable exten- 
sion of the luminosity range of X-ray starbursts previ- 
ously reported, but is consistent with the an extrap- 
olation of the X-ray-infrared relation for starbursts 
into the realm of hyperluminous infrared galaxies. 

7.3 Reliability of different luminosities 

Unless we happen to have 70 //m detections, estimates 
of starburst luminosities, and of infrared bolometric 
luminosities from Spitzer data arc uncertain by a fac- 
tor of ~ 2 (Rowan- Robinson et al 2005) . Because we 
are sampling the AGN dust torus template at its peak 
wavelengths, and often have several detected bands, 
the estimates of dust torus luminosity are more accu- 
rate, with an uncertainty ~ ± 0.1 dex. The starburst 
and dust torus templates are so different at 3-24 fim 
that there is little aliasing between them even with 
our minimum two detected infrared bands. 

What is the best way to estimate the total (X- 
ray to 3 fmi) bolometric luminosity of AGN ? For 



Type 1 objects the best estimate is undoubtedly from 
the extinction corrected optical luminosity, since this 
involves only a short extrapolation to the UV peak 
of the spectral energy distribution. We have used the 
optical- UV bolometric correction of Rowan-Robinson 
et al (2008), to give Lt^ ~ 2.0Lopt, with an uncer- 
tainty of ±0.1 dex. 

For Type 2 objects, or for Seyferts where the 
AGN is outshone in the optical by the parent galaxy, 
the absorption-corrected hard X-ray luminosity gives 
a less reliable estimate. We find a mean X-ray bolo- 
metric correction of 27, but with a big range, consis- 
tent with the 4-100 range found by Vasudcvan and 
Fabian (2007). We estimate the rms uncertainty in 
the X-ray bolometric correction to be 0.4 dex. 

For both Type 1 and Type 2 AGN we can use 
the luminosity in the AGN dust torus, assuming an 
average covering factor of 40%. The range of covering 
factors is 0.01-1, and the rms uncertainty in the esti- 
mate Lbh ~ 2.5. Ltor is 0.26 dex (Rowan-Robinson et 
al 2008). 

Finally it is worth remarking that our modelling 
of infrared SEDs will be enormously improved when 
Herschel data become available. 
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